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Abstract 

The single-electron energy and static charge-lattice deformation correlations have been 
calculated for the first excited state of a two-site Holstein model within perturbative ex- 
pansions using different standard phonon bases obtained through Lang-Firsov (LF) trans- 
formation, LF with squeezed phonon states, modified LF, modified LF transformation with 
squeezed phonon states, and also within weak-coupling perturbation approach. Compar- 
isons of the convergence of the perturbative expansions for different phonon bases reveal 
that modified LF approach works much better than other approaches for major range of the 
coupling strength. 
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I. INTRODUCTION 



The simplest model to study the nature and properties of polarons as a function of 
electron-phonon (e-ph) interaction is the Holstein model [|T| which consists of tight binding 
electrons coupled through a site diagonal interaction term to dispersionless phonons. The 
ground state properties of this model has been extensively studied during recent years. In 
the antiadiabatic limit the confinement of the lattice deformation around the charge carrier 
for large e-ph interaction gives rise to small polarons whose nature and dynamics is generally 
studied using the Lang-Firsov (LF) method based on the canonical LF transformation 0. 
For weak and intermediate e-ph coupling the LF method is not appropriate. For that region 
the importance and superiority of the modified LF (MLF) and the MLF with squeezed 
transformation (MLFS) over the LF method have been pointed out in previous works 
The results of exact diagonalization studies of a two-site Holstein model indicate the failure 
of the standard (zero phonon averaging) classical LF approach even in the strong coupling 
antiadiabatic limit and the authors of Ref . asserted that perturbation approach within the 
LF scheme is not meaningful. Marsiglio [0 studied the Holstein model in one dimension with 
one-electron up to 16 site lattices using numerical diagonalization technique and concluded 
that neither the Migdal approximation [|] nor the usual small-polaron approximation is 
in quantitative agreement with the exact results for intermediate coupling strength. So at 
present the Holstein model cannot be described by any single conventional analytical method 
for the entire range of the coupling strength either in the adiabatic or in antiadiabatic limit. 
In a recent work 0] we addressed this problem and considering a two-site system we developed 
a perturbation expansion using MLF phonon basis and the results obtained thereby are 
in good agreement with the exact numerical results for the entire range of the coupling 
strength for an intermediate value of hopping (t ~ uq). Subsequently, we presented a detailed 
comparison for the convergence of perturbation corrections to the ground state energy and 
wave function for the same system |]TU[ using different phonon bases obtained through the LF, 
modified LF (MLF) and modified LF with squeezing transformations (MLFS). The results 
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showed that for weak and intermediate couphng the pertubation corrections within the MLF 
and MLFS methods are much smaller and the convergence of the perturbation expansion 
is much better compared to the LF method. For strong coupling all the methods become 
equivalent and a good convergence in the perturbation expansion for the ground state is 
achieved. In this paper we have studied the convergence of the perturbation expansion 
for the energy and the correlation function for the first excited state of a two-site single 
electron system using different phonon bases obtained through the LF, MLF, MLFS, LF 
with squeezing (LFS) transformations, and also within weak-coupling expansion. 

In Sec. II we define the model Hamiltonian, describe different variational phonon bases 
states that we have considered and present the expressions for the energy, wave function and 
static correlation functions calculated within the perturbation method for the first excited 
state. In Sec. Ill we present the results obtained by different methods and discuss about the 
convergence of the perturbation series, hence the applicability of the methods in different 
regions of the e-ph coupling strength for different hopping parameters and conclusions. 

II. FORMALISM 

The two-site single-polaron Holstein Hamiltonian is 

i,a <y i,cr i 

where i —1 or 2, denotes the site. (c|^) is the annihilation (creation) operator for the 
electron with spin a at site i, and Uia- {=cl^Cio-) is the corresponding number operator, g 
denotes the on-site e-ph coupling strength, t is the usual hopping integral, bi and bj are the 
annihilation and creation operators, respectively, for phonons corresponding to interatomic 
vibrations at site i, and cuq is the phonon frequency. For the one-electron case spin index is 
redundant here, hence omitted in the following. 

The Hamiltonian (1) can be divided into two independent parts: one part describing sym- 
metric in-phase vibrations coupled to the total number of electrons and the other asymmetric 
out of phase vibrations coupled to the electronic degrees of freedom. The first part describes 
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just a shifted harmonic oscillator while the solution of the second part is a nontrivial problem 
[^,10 . The Hamiltonian for the second part is given by 



Hd = ^eni- t{c\c2 + c\ci) + ujQg+{ni - n2){d + d)) + UoSd , (2) 

i 

where = g/V^ and d = (bi — b2)/V^. 

For a perturbation method it is desirable to use a basis where the major part of the 
Hamiltonian becomes diagonal. The phonon bases chosen in the MLF or MLFS approach 
consist of variational displacement oscillators and can produce retardation effect even in 
the zeroth order of perturbation unlike the LF approach, hence these bases are better 
choices for perturbative calculation when the coupling strength is not very strong. We now 
use the MLF transformation where the lattice deformations produced by the electron are 
treated as variational parameters For the present system, 

Hd = e^//,e-^ , (3) 

where R = A(ni — n2){d^ — d), and A is a variational parameter related to the displacement 
of the d oscillator. The transformed Hamiltonian is then obtained as 

Hd = ujod^d + ^ epUi — t[ciC2 exp(2A((i^ — d)) 

i 

+ 4ci exp(-2A(rf"^ - d))] + uJo{g+ - A)(ni - n2){d + d^) , (4) 
where ep = e — uJo{2g+ — A) A. 



Now we will make a squeezing transformation [11] to Hamiltonian (4) 



H,d = e'HdC-^ (5) 

where 5* = a{didi — d\d\). This new phonon basis is squeezed with respect to the previous 
basis. The transformed Hamiltonian (5) takes the form 

Hsd = cc'o(i^c?[cosh^(2a) + sinh^(2a)] + uoq cosh(2a) smh.{2a){dd + rf^rf^) 
+ ^ epTii — t[c[c2 exp(2Ae((i^ — d)) + c^ci exp(— 2Ae((i^ — d))] 

i 

+ uJo{g-i- — X){ni — n2){d + d'^) exp(2a) + uoq sinh^(2Q;) (6) 



where Ae = A e . For the single polaron problem we choose the basis set 

|±,iV) = i=(et±4)|0)e|iV) (7) 

where |+) and |— ) are the bonding and antibonding electronic states and |A^) denotes the 
A^th excited oscillator state in the MLFS, MLF, LFS or LF bases depending on the method 
considered. It may be noted that the MLFS basis turns into the MLF basis when a = 0, 
and into the LFS basis if A = (7+. The MLF basis turns into the LF basis when X = g^, 
and it turns into the weak coupling expansion for A = 0. We consider the diagonal part 
of Hamiltonian (6) as the unperturbed Hamiltonian (-f^o) ^-iid the remaining part of the 
Hamiltonian, Hi = Hgd — Hq, as a perturbation [p|,p!0|. The unperturbed energy of the state 
|±, A^) is given by 



= (A^, ±\Ho\±, N) = cJo[sinh^(2«) + iV(sinh^(2a) + cosh^(2a))] 

TV ^2\ ,^i 



+ep =F te 



Ci 



(8) 



where te = t exp (— 2Ag). The general off-diagonal matrix elements of Hi are given in Refs. 



,|T0|. As noted in our previous work the unperturbed first excited state wave function 
should be built up as a linear combination of |+, 1) and |— , 0). The unperturbed energies of 
the states |+,1) and | — ,0) are (ep + cuQ — te(l— 4A^) +3sinh^(2a)) and (ep + tg + t^osiiih^(2a)), 
respectively within the MLFS approach, and the off-diagonal matrix element between these 
two states is (2Aete + '^0(5'+ — A)e^"). The unperturbed energies of these two states cross 
at an intermediate value of g+ for 2t > uq. So, following degenarate perturbation theory 
linear combinations of the states |-|-, 1) and | — ,0) are formed to obtain two new elements 
of basis states so that Hgd becomes diagonal in the sub-space spanned by these two states. 
The unperturbed first excited state is given by 

1^^)= , / [a|-,0) + &|+,l)] (9) 

V + ^ ) 

The ratio (c) of the coefficients 6 to a and the unperturbed energy (a) of the first excited 
state may be found out from the relation 



a — Hii 



H 
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(10) 



Hi2 01 — H22 

where Hu, H22, H12 are the matrix elements of Hd in the subspace of | — ,0) and 
Eq.(lO) gives two roots of a, the lower value of a (say, ai) corresponds to the first excited 
state. 

Our chosen phonon basis within the MLFS method have two variational parameters A 
and a, proper choices of them will make the perturbative expansion convergent. In our 
previous works [^,10| we find that the minimization of the unperturbed ground state energy 
of the system yields phonon bases (within the MLF and MLFS) for which the perturbative 
expansion for the ground state shows satisfactory convergence. The corresponding values of 
A and a are given by A = ^^_^^°^+_4^ and a = |sinh~^[2A((7+ — A)]. For the study of the first 
excited state we will use the same basis as we have used for the ground state. The first-order 
correction to the first excited state wave function is obtained as, 



1 



E 



Af=2,4,.. {Oil — J^-^N 
Wo 



\-,N) 



+ E 

N=3,5.. 



where We = W + VN cuoc(^+ - A) e^" 6^,2 + \ N{N - 1) ^ sinh(4a) Sn,2 



N) 



(11) 



Wo = W + y/NuJo{g+ - A)e2"5^,i + ^N{N - 1) f c sinh(4a) 5 

(2Ae)^ 



N,3, 



and W = te 



:i + 2AeC-|£). 



Second-order correction to the first excited state energy is given by. 



(2) 



1 + C2 



E 



+ E 



\Wo\ 



^A- ("1 - ^- 7v) Af=3,5.. ("1 - ^+ A 



(12) 



Using our previous prescription |TD| successive higher order corrections to the wave function 
d'^/'i^'')) well as to the energy [Ex^^) are calculated. 

The static correlation functions (niUi)o and (niM2)o, where u\ and U2 are the lattice de- 
formations at sites 1 and 2, respectively, produced by an electron at site 1, are the standard 
measure of polaronic character, and indicate the strength of polaron induced lattice defor- 
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mations and their spread. The operators involving the correlation functions may be written 
as 



Hi 




(13) 



The final form of the correlation functions are obtained as 



{niUi)o 



1 



{9+ + A) + 



(14) 



(^1^*2)0 



2 



2 



1 



{9+ - A) 



Ngi 
Ngi 



where A^'gi = {ipi\ipi), Ai = {ipilrii^d + S)\ipi) and \ipi) is the perturbed wave function of 
the first excited state. 



The perturbation corrections to the energy and wave function of the first excited state 
are estimated up to the sixth and fifth-orders, respectively, within the LF, LFS, MLF, MLFS 
methods, and weak- coupling perturbative method considering 25 phonon states (which is 
sufficient for < 2.2) in the transformed phonon basis. 

In Fig. 1 the perturbation corrections to the energy of the first excited state are plotted 
as a function of for t/ujQ = 1.1. The perturbation corrections in energy are smaller within 
the MLF and MLFS than the LFS and LF methods. For the LF, LFS, and MLFS methods 
the convergence of the corrections is weaker in a range of g+ where third and fourth or fifth 
and sixth-order corrections are comparable. Within MLF approach energy perturbation 
corrections show satisfactory convergence in the entire region of (7+ and the energy of the 
first excited state, when computed considering up to the fifth or sixth-order corrections, 
becomes identical with the exact result of Ref. [^]. 

The shape of the correlation function (^1^2)0 is much more sensitive to the corrections 
to the wave function than {niUi)^ which indicates that convergence of {niU2)o is a clear 
signature for convergence of the wave function. In Fig. 2 we plot the correlation function 
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(niM2)o obtained by considering up to the different orders of perturbation corrections to the 
wave function against g^. The LF method shows a bad convergence for low values of (7+, 
while a good convergence beyond (7+ = 1.0. The convergence within the LFS method is not 
satisfactory over a wide region of (7+. The MLF method shows very good convergence for 
low values of (< 1.0) as well as for high values of g^ (> 1.3). Convergence within the 
MLFS is excellent for low values of (7+ , but then for a wide region of intermediate values of 
(7+ the convergence is not satisfactory. When compared with exact results of {niU2)o (taken 
from the Ref. 0), it is found that the MLF results up to the fifth order perturbation are 
identical with the exact results except in the region 0.9 < g+ < 1.3 where a slight departure 
in values from the exact results is seen. 

In Fig. 3 we plot the perturbation corrections of different orders to the first excited state 
energy and the the correlation function {niU2)o evaluated by considering up to different 
orders of perturbation corrections to the wave function against (7+ within the weak-coupling 
perturbative expansion for t/ujQ = 1.1. The odd order corrections to the energy within 
the weak-coupling expansion are zero. The energy corrections of any (even) order increases 
monotonically with g^, since the entire e-ph interaction is treated as a perturbation within 
the weak-coupling scheme. The correlation function within the weak-coupling perturbation 
procedure shows good convergence for low values of (7+, as expected. However even in this 
region, convergence within the MLF and MLFS methods is found to be slightly better than 
that obtained within the weak-coupling scheme. 

In Fig. 4 we show the energy corrections and correlation function, calculated by consid- 
ering different order of corrections, within the MLF method for t/ujQ = 0.6. The figure shows 
that the perturbation corrections of successive orders are really very small and shows perfect 
convergence for the entire region of g+. The correlation function calculated considering up 
to second order correction to the wave function would reproduce almost the exact result in 
this limit. For lower values of t obviously the convergence within the MLF scheme would 
be much better. We have also studied the perturbation corrections to the energy and the 
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correlation function {niU2)o for the first excited state in the adiabatic region of hopping 
{t/uo — 2.1) within MLF approach and found that energy convergence is quite good for 
the entire range of whereas convergence of corrections to the wave function is fairly well 
except for intermediate range of coupling. 

In conclusion, our study on the first excited state of a two-site single electron Holstein 
model shows that the perturbation method based on the MLF variational phonon basis is 
better than the other methods (LF, LFS, MLFS) if the entire range of g+ is considered. 
The MLF method could yield exact results for the entire range of g+ in the antiadiabatic 
cases {t/ujo < 0.6) and for a major region of g+ (covering both low and high values of 
g+) for intermediate hopping {t/ujQ = 1.1). For small values of g+ both the MLFS and MLF 
methods show excellent convergence, convergence is found to be even better than that within 
the weak-coupling perturbation method. Our previous study on the ground state and the 
present study on the first excited state thus establish that the perturbation method based 
on the MLF phonon basis could satisfactorily describe the ground state as well as the first 
excited state of the two-site system for major region of the e-ph coupling strength. 
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Figure captions : 

FIG. 1. Variation of the perturbation corrections e[^^ to the first excited state energy as 

a function of the couphng strength (g^) for t/ujQ = 1.1 in (a) LP, (b) LFS, (c) MLF, and 
(d) MLFS methods. e["'^ is the nth order perturbation correction to the first excited state 
energy. 

FIG. 2. Plot of the correlation function (^11*2)0 calculated up to different order of pertur- 
bations in the wave function vs g+ for t/uJo = 1.1 within different methods (a) LP, (b) LPS, 
(c) MLF, and (d) MLFS. The labels (2), (3), .... denote the curves obtained by considering 
up to the second-order, third-order, .... corrections to the wave function, respectively. 

FIG. 3 Convergence within weak-coupling expansion for t/cuo — 1.1 : (a) the perturbation 
corrections to the first excited state energy and (b) correlation function {niU2)o calculated 
up to different order of perturbations in the wave function, as a function of 

FIG. 4 Convergence within MLF approach for t/cuo — 0.6 : (a) the perturbation correc- 
tions to the first excited state energy and (b) correlation function (7111*2)0 calculated up to 
different order of perturbations in the wave function, as a function of g+. 
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FIG. 3 
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Fig. 4 



